Separability and entanglement in finite dimer-type chains in general transverse fields 
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We determine the conditions under which general dimer-type spin chains with XY Z couplings 
of arbitrary range in a general transverse field will exhibit an exactly separable parity-breaking 
eigenstate. We also provide sufficient conditions which ensure that it will be a ground state. We 
then examine the exact side limits at separability of the entanglement between any two spins in 
a finite chain, showing that in the vicinity of separability, the system will loose all signatures of 
dimerization, with pairwise entanglement approaching infinite range and becoming independent of 
separation and interaction range. The possibility of a non-uniform exactly separable ground state 
induced by an alternating field is also shown. As illustration, we examine the behavior of the 
pairwise entanglement in a finite XY dimer chain under a uniform as well as alternating field. 
Related aspects of the magnetization are also discussed. 

PACS numbers: 75.10. Jm,03.65.Ud,03.67.Mn,64.70.Tg 



I. INTRODUCTION 

Quantum entanglement is an essential resource for 
quantum information science, allowing radically new 
forms of information transmission and processing^—. It 
has also aroused great interest in condensed matter and 
many-body physics^, providing a novel perspective for 
the analysis of strongly correlated systems. Fundamental 
properties of entanglement in quantum spin chains have 
been determined, especially in connection with critical 
phenomena in the thermodynamic limil£~— . The study 
of finite spin chains can also provide new insights into 
the most basic aspects of entanglement, and is presently 
also stimulated by the unprecedented level of control 
that can be reached in some recently developed quantum 
devices^, able to realize spin arrays with controllable 
Heisenbcrg interactions. 

A fundamental related question is the range the en- 
tanglement between individual spins can reach under the 
action of an applied magnetic field. At the standard crit- 
ical field of large anisotropic XY or XY Z chains, it re- 
mains finite and typically small (for instance, restricted 
to just first and second neighbors in a II? Ising chain in a 
transverse field 6 ). However, it can diverge at a different 
field: Anisotropic chains may also exhibit a factorizing 
field, where an exactly separable ground state (GS) be- 
comes possible, i.e., where the mean field GS becomes 
exact. This remarkable feature was first discovered in 
ID chains with first neighbor coupling s 10 ' 11 , and recently 
examined in detail in more general systems in a uniform 
fielder— . A general method for determining separability 
was in particular developed in refs i 17 ' 18 . In the immedi- 
ate vicinity of the factorizing field, the pairwise entan- 
glement in a finite chain can reach full rang o 14 ^ 16 . The 
transverse factorizing field in finite XY Z chains arises ac- 
tually at the crossing of opposite S^-parity levels^ 6 -, with 
separable parity breaking eigenstates emerging from the 
superposition of the entangled definite parity states^ -. 

The aim of this work is to examine the previous issues 
in finite dimer-type arrays, which have recently received 



much attentio n 17 ' 19 ' 21 ' 22 . We will consider arrays of ar- 
bitrary spins with XY Z couplings of arbitrary range in a 
transverse field, not necessarily uniform, and determine 
the separability conditions together with the entangle- 
ment side limits at separability, which will be shown to be 
independent of separation, coupling range and other de- 
tails such as the strength of the coupling between dimers. 
At these points all traces of dimerization will then be lost. 
We will also examine factorization under an alternating 
field, which can give rise to a separability curve with field 
dependent separable solutions and entanglement limits. 
Entanglement between spins unconnected by the inter- 
action can in this way exceed that between linked pairs 
in the vicinity of separability. These effects are specially 
noticeable for finite chains close to the XXZ limit. Re- 
lated aspects of the magnetization and the entanglement 
between one spin and the rest of the chain are also dis- 
cussed. 

Section II contains the general theoretical results, in- 
cluding the mean field+RPA interpretation of the separa- 
bility conditions, while sec. Ill the application to general 
dimer-type chains, including illustrative exact results for 
finite chains. The appendix contains the details of the 
exact calculation obtained through the Jordan Wigner 
transformation. Conclusions are finally drawn in IV. 



II. FORMALISM 



A. Transverse Factorizing Fields 



We first consider the general Hamiltonian 



S i S j + V V S i S j 



'), (1) 



which describes an array of n spins Si not necessarily 
equal, interacting through XYZ-type couplings of arbi- 
trary range in a general transverse applied field 6\ It 
satisfies [H, P z ] = 0, where P z = exp[i7r Y^i=i ( s f + ,s «)] 



2 



denotes the global S z parity or phase-flip (here is 
the spin value at site i). Denoting with |0j) the local 
state with its spin fully aligned along the — z direction 
(sf|0j) = — Sj|0j}), this Hamiltonian will exhibit a fully 
separable parity breaking eigenstate of the form 



|6) = \e, 



=1 exp[zt 



»?]|0i 



(2) 



i.e., a state with its spins fully aligned along local axes 
forming angles 9i with the z axis, if (and only if) the 
conditions 



v y = v x cos Qi cos Qj + yl z sm 9% sm $j 



(3) 



6 l sin^ = 



> s j v* J cos 0i sin #j — u*- 7 sin #j cos #j ) ,(4) 



\ " 



are satisfied 2 ^.. They can be obtained replacing in 
H by the rotated operators e~ iSi e iSi and solving 
Hq\Q) = E \O), where Hq = e-^^Wffe'Si^? and 
|0> = ®i|0i). Eqs. ©-© actually hold for general local 
rotations e l( "' Si |0j) since the latter can also be cast in 
the form ([2]) through complex angles 9i and a suitable 
normalization facto r 20 i 23 . Note also that for a spin 1/2 
array Eq. © is in fact the most general separable state. 
The energy Eq becomes 

n 

Eq = - s i W cos $i + | s j ( v x sm Oi sm &j 



i=l 



cos Oi cosOj)] 



(5) 



If l^'l < v l J and 6>j G (0,tt) V i, |0), as well as 
its degenerate partner state 

\-e) = p z \G) = \-e 1 ,...-e n ), 

will be ground states of H when Eqs. ©-© arc 
fulfilled 20 . Of course, they can be GS also in other 



case; 



,10.11 



by suitably adjusting the relative signs of the 
9i (see sec. MI A[) . Eqs. © are in fact the stationary 
conditions for the energy © at fixed b l , v 1 ^ , representing 
the mean field equations which ensure stability of | ± 0) 
against one-spin excitations. 

Eqs. ((3]) warrant that |±0) will be exact eigenstates by 
canceling the residual matrix elements linking |±0) with 
two-spin excitations, and have a clear meaning within 
the random phase approximation^^- (RPA): If satisfied, 
Vi,j the RPA vacuum will coincide with the mean field 
state. More explicitly, the zero temperature RPA matrix 
(whose eigenvalues are the RPA energies) is 



■Hrpa = 



Br 



A 

-B~ 



B- 



-A 



■W 



V~ J smt/i sinf 



where ±^Ai are the eigenvalues of the local mean field 
Hamiltonian b l sf — V - 
ments associated with the dispersion {s^s~) and creation 



( s j)eSi, and Bf- the ele- 



(sfs^) of spin excitations respectively. Eq. (J3j) is then 
equivalent to the condition B~ =0, implying no RPA 
corrections to the mean field vacuum. 

From Eq. © it is seen that a uniform eigenstate with 
Oi = 9 Vi becomes feasible if the anisotropy ratio 



_ v 
X a 



cos 2 9 . 



(G) 



is constant for all pairs and satisfies x > (if x > 1 
(complex 9), a global rotation of ir/2 around the z axis 
will lead to x — )• 1/x and 9 real). Eq. © leads then to 



(7) 



if x G [0,1) (the opposite sign for all b % is obviously also 
feasible) and to b % arbitrary if % = 1 (XXZ or XX 
case2£, where |0) is an exact eigenstate V V). Any spin 
array with couplings satisfying Eq. © will then exhibit a 
uniform separable degenerate eigenstate |0) = \9, . . . , 6) 
if the fields b % at sites i are tuned at the values ©. It 
will be a GS when \vt j \ < v l J Vi,j. 



B. Entanglement at factorizing fields 

In a finite array the exact GS of H will not be in general 
exactly degenerate away from the factorizing point, and 
will have a definite S z -parity. The correct side limits at 
the factorizing field are then provided by the normalized 
definite parity states 



l© ± ) = 



|e)±|-e) 

V2(l±(-0|0»' 



(8) 



where (— 0|0) = T\™ = i cos 2si 9i is the overlap between 
the degenerate separable eigenstates. The states © sat- 
isfy P z |0 ± ) = ±|0 ± ) and are obviously also exact eigen- 
states when Eqs. ©-© are fulfilled. 

These states are entangled, with Schmidt ranki 2 for 
any bipartition (A, A) of the whole system^ (here A de- 
notes a subset of spins and A the complementary sub- 
set). Moreover, the reduced state of any subsystem of 
two or more spins can be effectively considered as a two- 
qubit mixed state with respect to any bipartition^S. The 
entanglement between any two subsystems can then be 
measured through the concurrence, a measure of entan- 
glement originally introduced for two qubit systems^ 
(where it can be exactly computed, see sec. MI D[) . and 
later extended to mixed states of general bipartite sys- 
tems through the convex roof extension of the generalized 
pure state expressio n 2 ^ 29 . The concurrence between any 
two spins i,j in the states |0 ± ) can be shown to be^ - 



cos* 



4s.,- 



i±(-©|0) 



(9) 
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where (— Qjj\Qjj) — Jlfc/i j cos4 ' £ denotes the comple- 
mentary overlap. It will be appreciable for sufficiently 
small angles 9 k if 6i,0j ^ 0. 

On the other hand, the entanglement between one spin 
and the rest of the chain in the states ((8]) can be mea- 
sured through the entropy Si = — Tr[p^ logpi], where 
Pi = T?i |0 ± )(G ± | is the reduced density matrix of the 
spin at site i, or alternatively, through the pure state 



concurrent 



,2S 



C, = y/2(l - Trpf). The latter provides 
an upper bound to the sum of all pairwise concurrences 
dj stemming from site Cf > G%. It fully 

determines Si when Si = 
in the states (El) reads 20 



1/2 (sec. IIII D|) . Its expression 



C 



(i - cos 4 ** 00(1 - c ° s4sfc °k) 



1±(-9|9) 



(10) 
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FIG. 1: (Color online) Schematic plot of a) a system de- 
scribed by the couplings {TTJ, representing two interacting 
cyclic chains and b) the dimer chain of Eq. (|12[) . a particular 
case of a). Numbers indicate the notation used in figs. [2H31 



with Cf w \/l — cos 4si 9i if the overlap is neglected. The 
entanglement between L and n — L spins, as well as be- 
tween any two sets of spins in the states (EJ can also be 
exactly calculated^. 

When | ±6) are GS, Eqs. ©-(fTUl) represent the actual 
side limits of the GS concurrences CV,- and Ci at the fac- 
torizing point, where a transition |0 + ) — > |0~) will take 
place as the field increase s 16 ' 20 . The entanglement be- 
tween two spins will then reach full range in its vicinity, 
provided and fl 4 ^ 0, V i,j (Eq. ©). 

When 6i = tt/2 Vi (corresponding for v z = to the 
Ising case v y — according to Eq. (EJ)), | ± 9i) are or- 
thogonal and Cf~ — 1 while Cy = V i, j. The previous 
effect becomes significant in the opposite limit of small 
6i (systems close to the XX Z limit). We also remark 
that the uniform mixture of both definite parity states, 
pa = i(|e + )(6 + | + |e-}(6-|), is also entangled and 



leads to attenuated concurrence; 



: 20 Cfj 



C 



c? = c. 



-&m 



<-e|e) 
ij i+(-e|e> ' 



i i+<-e|e>- 



III. APPLICATION TO DIMER-TYPE CHAINS 

Let us now consider a pair of uniform interacting chains 
of the same size m and spins s a , not necessarily equal. 
We can embed this system in a single non-uniform chain 
of even size n — 2m assigning odd (even) sites to the first 
(second) chain, as schematically depicted in Fig. [T] (left), 
such that a = o, e. We may then consider 



i) 



(11) 



where at — o, e indicates the parity of the site, such that 
v"(l) represents the interchain (internal) couplings for I 
odd (even). Accordingly, v°(l) — v e (—l) for I odd and 
v a {l) = v a (-l) for I even (but v°(l) f v e (l) in general). 
In the cyclic case v a (— I) — v a (n — I) V a, I. 

An important example of this type is that of a dimer 
chain with just nearest neighbor couplings (Fig.[Tl right), 



where <(/) = v^S L±1 : 



H,, 



n/2 



i=l 



E/ 

li—x,y,z 



V u S 2i-l S 2i +V u S 2i S 2i+ll 



(12) 



Here v e ^ can be considered as the (weak) couplings be- 
tween dimers and w° the (strong) internal couplings, the 
system becoming dimerized (i.e., an array of independent 
spin pairs) for — > (see also sec. Ill C). 

A different example of (fTTj) . which nonetheless will ex- 
hibit factorization properties similar to those of Eq. ([12")) 
(see below), is a pair of arrays with no internal cou- 
plings interacting through a constant full range coupling: 
v°(l) = 2v^/n V I odd and otherwise, such that 

H p = b°Sl + VSt-± J2 v»SZS£, (13) 



where S% e = ^ odd sf are the total spin components of 
each array and we have assumed a constant field is ap- 
plied to each of them. This system is obviously equivalent 
to an interacting pair of spins S — \ns and S e = \ns e 
if restricted to the maximum spin multiplet. As in the 
Lipkin model^ 2 -, the 1/n scaling ensures here a bounded 
intensive energy (H) /n for n — > oo and fixed . 



A. Uniform separable eigenstates 

In the general case (TTTj) with cyclic boundary condi- 
tions, a separable eigenstate with a common angle 9i = 9 
\/i is then feasible if Eq. (E|) holds for any connected pair, 

i.e., x = V J" (i)- v Z " (i) e [0' 1) an( l constant VZ, and there is 



a uniform field b 1 — b ai in each subchain given by 



E 



(14) 
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where «™ - J2i even <(0 " <(0, = odd <(0 " 
v^(l) — v eo , with 6° = 6 e if s D = s e . Such uniform 
eigenstate is also feasible for a similar open chain pro- 
vided a non- uniform field, as determined by Eq. ([7]). is 
applied. For short range couplings this will imply just 
boundary corrections. The ensuing states | ± 0) will be 
GSif \v%(l)\ <<(/) Vcr,Z. 

The definite parity states ([8]) will then lead to a finite 
concurrence © for any spin pair, which will depend on 
the parity of the sites but not on their separation: The 
odd-odd (C^,), even-even {Cf e ) and odd-even (Cf e ) con- 
currences will be given, according to Eq. ([5]), by 



{i-x 2s °)x 



,S-2s„ 



ct 



VcEcE, (is) 



where S = ^n(s a + s e ) is the total spin. The range of 
the entanglement between two spins will then increase 
as the factorizing fields (IT41 are approached in each sub- 
chain, becoming independent of the coupling range and 
separation. If s a = s e , obviously C 00 — C ee — C oe . 

C^ a will be appreciable for sufficiently small XY 
anisotropies: If x ~ 1 — then x S ~ e ~ S for small S/S. 
In fact, for S -> (XX limit) C+ CT -> but C~ a -> 2s ff /5 
(i.e., 2/n for s Q = s e , which is the maximum attain- 
able value when all pairs are equally entangled^), as 
|9 + ) — s- |0) but |0~) approaches the entangled W-type 
stated ex J2i s t\®)- I n tne opposite limit (Ising case 
X = 0), | ± 0) become orthogonal and C aa = 0. 

In the dimer chain (|12[) . the uniform separable eigen- 
state becomes then feasible if there is a common 
anisotropy x — V Z _ v l € [0, 1) for a = o, e, and the fields 



are chosen as (Eq. (Q~4 



b° = Jxv° 



b°s /s e 



(16) 



where v oe — 



In an open chain we 



should just add the border corrections b 1 — \b° , b n — }-b e 
according to Eq. ©. Thus, in the ferromagnetic-type 
case 1 < v% for a — o, e, its GS will become uniform 
at the factorizing fields (|16l) . regardless of the ratio v^/v° 
(as long as it is non-zero), loosing there all signatures 
of a dimerized structure and leading to the full range 
concurrences ([TBI as side-limits. 

Let us also remark that for the nearest neighbor cou- 
plings of (fT2"l) . the antifcrromagnetic case v% < Vct 
can be brought back to the previous case by means of 
local rotations of angle 7r around the z axis at even 
sites (implying sf — > (— l) l+1 sf and hence v° — > —v° 
V a for « = x.y). A uniform separable eigenstate 
\/ A / A . . .) = \06 . . .) in the rotated system corresponds 
then to an alternating solution 9i — (— l) l+1 9 (Neel-type 
state |/ , \^\ . . .) = \6, -9, 9, . . .)) in the original sys- 
tem. Note that this holds for arbitrary spins s a (equal 
or distinct). Separability (but not entanglement) in the 
s o = s e — 1/2 dimer chain was discussed ir*i£, with the 
correct treatment for general antiferromagnetic couplings 
discussed in detail in——*-". 



For even m = n/2 (to avoid frustration effects^), the 
mixed case v° > 0, v% < (or viceversa) can also be recast 
as a ferromagnetic case v% > Vct by means of local rota- 
tions of 7r around the z axis in even sites of both subchains 



( S 2i+fc 

ing v e 



(-1) 



"2i+k 



for k — —1, and ji = x,y, imply- 



-Vp). The uniform solution corresponds here 
to f?2 S :+/c = (— l) l+1 9 for k = —1,0 in the original mixed 



system, i.c 



18,19 



...) = \t 



Hence, for even n/2 we may always assume v° > Vct in 

In the system (|13p . the same uniform separable eigen- 
state becomes feasible if x = ^-1* e ano - ^ e 
fields are set at the values (fl6|) . with w oe = ^ — v z (Eq. 
([HI) )- We may again assume v x > since its sign can be 
changed replacing S£ — > —5^ for /i = x, y. This system 
will exhibit just three different GS pairwise concurrences 
(C 00 , C ee , C oe ) for any 6°, 6 e , which will approach the 
same limits IT5l) at the fields (1161). 



B. Alternating separable eigenstates 

In the case of two interacting subchains with no inter- 
nal couplings, like Eqs. (jT2"j) or (fT5)l . we may also consider 
the possibility of different and controllable uniform angles 
9 Q , 9 e (with \9 \ 7^ |0 e |) in each subchain, i.e., 

|0) = \9 9 e 9 9 e ...), 

through an alternating field b° ^ b e . For simplicity we 
will consider XY couplings (v l J = 0). According to Eqs. 
©"©J such a solution is feasible if for u = o, e and I 
odd, 

cos0 o cos0 e , b°=v oe ^fs S} (17) 



tan 0^ 



where w oe = ^ z odd v%(l) and a ^ a (i.e. w oe = u° + v^; 
in (JUJ) and w oe = w x in (fT5|)). This implies 

,.oe\2. 



b°b e = x(v 0e YSoSe, 



cos 



r„ i Oct — 
1 + 02 



(18) 
(19) 



Hence, for fields b°, b e satisfying Eq. (|T8|) we obtain a 
separable eigenstate with alternating angles 9 Ql 9 e deter- 
mined by Eq. (|19p . Since one of the fields is now free, 
we may adjust in such system the individual angles and 
thus the internal (Cf a , Cf e ) and interchain (Cf e = Cf a ) 
pairwise concurrences at separability, given now by 



(1 



-\ / 2 s £j \ •v/tS'o- 2s CT -w'S'f 
Act /Act Act 



l±X^xf s 



C± = y/C&C&(?D) 



where Xa = cos 2 f9 CT and 5^ = ns a /2. If |6°| > |o e | and 
s = s e , Cf a < Cf e < Cf e , despite the absence of even- 
even direct coupling (v ee = v°° = 0). For b° = b e we 
recover Eqs. <[T5]HT6l). The values of C± a , depend now on 
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the ratio rj = b°/b e (b a = y/rjx, b e — \fxpq when Eq. (fT8)) 
holds). For r) >■ 1, cos# — > 1 but cos# e — > x ; implying 
that in this limit remains finite at the factorizing 
field, while C^ D and C^, vanish. Note also that 6* CT is a 
decreasing function of b a . 

In the ferromagnetic case v% > 0, 9 and # e have both 
the same sign. For antiferromagnetic couplings v a x < V 
a in the dimer chain (fT2|) . we would have instead O > 
and 9 e < (or vice- versa), whereas in the mixed case 
v%v% < 0, 1 8) = \0 o e ,-9 o ,-6 e ,9 o 9 e ,...), as previously 
discussed. Border corrections b 1 = hb°, b n = i& e would 
also apply in an open dimer chain. 

The concurrence between one-spin and the rest of the 
chain C ai = Q, will be given at separability by (Eq. (|T0|) ) 



array of maximally entangled pairs) at zero field, since 
in this case Cf 2 = 1 and all eigenstates ) are Bell 
states. However, at finite fields, Cf 2 = 1 only if b± = 
(b e = =f6°), in which case just half of the eigenstates 
remain maximally entangled. For b° ,e > 0, maximum 
entanglement (C12 = 1) for — > will then arise just for 
an homogeneous field b + = b lower than the factorizing 
field, i.e., when the pair GS is antiparallel (|^^)). 
Let us finally notice that for v = ±, 

= (ni<in) = + v -» w 

implying (sf)+ = (s|)+ but (sf)_ = -(sf)-> Le - oppo- 
site magnetizations for negative ^-parity (see below). 



D. Results 



C. Spin 1/2 pair 

We may explicitly verify the previous expressions 
(valid for general n) in the "two qubit" case (s a — s e = 
1/2, n = 2), which also represents the — > limit in the 
spin 1/2 dimer chain (fT2"|) . Setting v± — \{v x ± v y ) > 
and 6± = \{b°±b e ), with 6° = 6 1 , b e = b 2 , the eigenstates 
and energy levels of Hamiltonian (JlJ become in this case 

= a" Itl) ± aillt) , B£ = K T ^ + ^(22) 
= a+|tt) ± a+U) , E+ = -\v z T 

where (a±) 2 = |(1 ± ) and the superscript 

i/ = ± indicates the S^-parity. The GS corresponds to 
|*+) or |*|), with E+ and £| crossing precisely when 
the factorizing conditions hold. At this point 

|*+) become the states @. In particular, for an ho- 
mogeneous field (6_ = 0, b + = b), ET = E+ when 
b = \yjx( v x — v z ) (Eq. (UJ)), whereas for w z = they 
cross when b°b e — \x v x (Eq- (|T5|> ). It is then explic- 
itly verified that the states (J5]) are the true side lim- 
its at the crossing point, with separability arising just 
from the crossing of these two states. Factorization cor- 
responds then to the quantum critical point of the spin 
1/2 pair. It should be also noticed that |*+) can here 
be always written as projected states ([2]) using suitable 

angles (tan 2 h6\ = °7°: , tan 2 h9 2 — a - a + ), The con- 
currence C12 = y / 2(T = Trp 2 ) in the states |*^) reads 

C± = 2|a±e£| = KI/^+4, (24) 

and coincides with both general results ([^-(ITUl) for the 
present case {Cf 2 = |sin0i sin#2|/(l ± cos#i cos #2))- 

In the spin 1/2 dimer chain (|12l) . Eq. (j2~4")l represents 
the limit of the concurrence C^i—i^i for v e ^ — ¥ 0. This 
implies that its GS will become fully dimerized (i.e., an 



Figs -HHS1 depict illustrative results for the GS pairwisc 
concurrence Cij in a finite spin 1/2 dimer chain described 
by Eq. (fl"2"|) with cyclic XY couplings (<y£ = 0). We have 
set = cwZ, with v° = x v ° x f° r = o,e. Full exact 
results for finite n can in this case be obtained through 
the Jordan- Wigner transformation (see Appendix). In 
this system Cij — Cij-i+i (Czj-i+i) for i odd (even). 

The reduced density matrix for a pair of spins i,j will 
commute with the pair parity e l7 '^ Si+s i~ 1 \ being then of 
the form ((. . .) denotes here the GS average) 

Pa = \ + <*?K + + 4 E «4>f s " ( 26 ) 

The GS pairwise concurrence Cij can then be evaluated 
as2£ 2A max — Tr i?, with A max the greatest eigenvalue of 

the matrix R = 4^/ sf s V j pij sf pij , and reads 

Cij - Max[C+,CT,0], (27) 
C± = 2[|< S f S | T 4^)1 - y/Tl T - i( S f T S , 2 ) 2 ] , 

being parallel^ (i.e. , as in the states |ft> + ) if Cf- > 
and antiparallel (|tl> - lit)) if > (just one of 
can be positive) . The entanglement of formation of the 
pair can then be obtained as^£ S{j — — J2v=±Pv^°S2Pv, 

where p± = i(l ± Jl - C?.). Gy = % = (1) for a 

separable (maximally entangled) pair. 

The case of a uniform field b° = b e = b is depicted in 
Fig. [5] Here Cij = C2.J+1 for j odd. At b — and for a = 
1 (uniform chain), there is entanglement just between 
first neighbors (C12 = C23 > 0). For a small anisotropy 
X = 0.9, as soon as the ratio a decreases below 1 the 
concurrence between weakly coupled pairs (C23) rapidly 
decreases (top panel), vanishing here already for a < 
0.74, whereas Cyi rapidly increases, practically reaching 
saturation for a = 0.25 (center panel). Hence, at zero 
field approximate dimerization is achieved already for low 
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b°/v x ° 



FIG. 2: (Color online) Concurrences between spins i,j vs. 
magnetic field in a spin 1/2 XY dimer chain (Eq. p2p ) for 
two values of a = vt/vZ and of the anisotropy \ = v y/ v x- 
The field is here uniform, with n — 20 spins. All C\j approach 
the same side limits (|15f) (which are independent of a) at the 
factorizing field (|28|l (red bar), as seen in the insets (blow 
up of main plot), changing from antiparallel to parallel and 
exhibiting there the same finite step. The red dot at b — b s 
indicates the concurrence Cfj in the mixture of both definite 
parity ground states. 



a, the system becoming essentially an array of maximally 
entangled pairs in the antiparallel states I^X)- 

The previous picture remains valid for weak finite 
fields. As seen in the top and central panels, increas- 
ing the uniform field destroys dimerization in a stepwise 
manner, the GS remaining almost unchanged until the 
first step, occurring at b* sa hVx v xQ ~ a ) f° r a n °t 
close to 1. These steps, clearly visible in small chains 
with low anisotropy, reflect the n/2 GS SVparity tran- 



sitions (crossings between the lowest levels of opposite 
parit y 16 ' 19 , which are close but not degenerate) taking 
place as the field increases when x € (0, 1], as in the ho- 
mogeneous XY chain. At the same time, the concurrence 
range increases as the last step is approached. The latter 
occurs precisely at the uniform factorizing field (Eq. ([7])) 



(28) 



where the dimer structure is completely lost and entan- 
glement reaches full range: All pairs become equally en- 
tangled irrespective of separation or location, with Cy 
reaching the side limits (j 15[) V i ^ j (lim 



which are independent of a and hence the same in top 
and central panels. At this field all Cy exhibit the same 
finite discontinuity, changing from antiparallel (b < b s ) 
to parallel (b > b s ). For a 0, b* and b s (first and last 
steps) merge at the two-qubit factorizing field \^[x v °x- 

For stronger fields b > b s we obtain a weak parallel 
concurrence, which for first and second neighbors persists 
for arbitrarily strong fields and can be described pertur- 
batively. First (second) neighbors concurrences are first 
(second) order in v%/b and given, up to 0(v°/b) 2 , by 



C 12 « l^|-3(«^) 2 . ^3 



if — ) 2 



i,i+2 



bb° 



±|_l(^)2 (1+a 2 ); 



where v° ± = \{v° x ± v° y ) = \v%(\ ± X ) and b = + 6 e ). 
Note that a threshold value of a is required for a positive 
second neighbor concurrence for strong fields. 

For higher anisotropics (lower x), the behavior be- 
comes similar to that of larger systems. The GS parity 
transitions become less noticeable, as seen for \ = 0.5 
in the bottom panel, and the pairwise concurrence side- 
limits at the factorizing field are smaller. Nonetheless, 
the increase of the concurrence range in its vicinity re- 
mains clearly appreciable. Let us remark that for small 
separations \i—j\, the results for for n — 20 at \ = 0.5 
are already very close to those for n — > oo. We should 
also mention that as x decreases, lower ratios a are re- 
quired to achieve approximate dimerization at low fields 
(at x = 0.5 and b = 0, C23 vanishes only for a < 0.58). 

Fig. 3 depicts the typical behavior for small 
anisotropics when different fields are applied at even at 
odd sites, with a fixed ratio rj = b°/b e = 3. The factoriz- 
ing value for b° (Eq. (|18[) ) is here 



(29) 



At b° there are now three different limits for the concur- 
rences at each side, Cg , Cf e and C^ e , which represent the 
common side limits of Ci ; 2j+i, ^2,2^+2 and C\^j, ^2,2^+1 
Vj and are given by Eq. ([20)) . They satisfy here 



Cte/Cto = (x + ?l)/(x + ri- 1 )>l 

for 77 > 1, implying Cf e > Cf e . In particular, C24 (which 
approaches Cf e for b° — >• (b")^) clearly exceeds in the 
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FIG. 3: (Color online) Same details as the top and center 
panel of fig. [5] for the case of different fields at even and odd 
sites, with a fixed ratio b° /b e = 3. Now odd-even, odd-odd 
and even-even concurrences approach different common side 
limits at the factorizing field (|29[1 . with C24 becoming the 
greatest concurrence in its vicinity. 



vicinity of b° both first neighbor concurrences C%2 and 
C23, despite the absence of second neighbor couplings. 

It is also seen that C\i is no longer nearly constant up 
to the first parity transition, which occurs now at b°* w 
hy/VXVxi^ ~ a )- This effect can be easily understood 
with the two qubit concurrence (|24|) : At low a, C12 is 
essentially described in the first region by the two qubit 
expression flM|) for C~ , which for &_ = \b°{\ - 1/rj) ^ 0, 
is no longer constant and decreases as b° increases. 

In order to highlight the universality of the limits at 
the factorizing field, wc depict in fig. 2] the pairwise con- 
currences in the system (|13l) for s Q = s e = 1/2, where 
the exact solution can be obtained through direct diag- 
onalization in the S Q = S e — n/2 representation. In the 
XY case, v oe = v x and the odd factorizing field at fixed 
ratio b°/b e = rj is 6° = \v x ^/rfx- There are now just 
three different concurrences at all fields, C 00 = Ci^j+i, 
C ee = C 2 ,2j+2 and C oe = C\^j = C eo ( j-independent) , 
which approach the same limits of figs. (Eqs. (|T5|). 
(|20|l) at the factorizing field, since the latter depend 
solely on \ an d the field ratio r\. They are here compara- 
ble to the values away from the factorizing field, since the 
monogamy bound on Y^j^i Ch entails C aa < — 0(1/ n) in 
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FIG. 4: (Color online) Concurrences between spins i,j vs. 
magnetic field in a system with constant full range couplings 
between even and odd sites, described by Hamiltonian (|13|) . 
We have set again \ — 0.9, n = 20 and a uniform (alter- 
nating) field in the top (bottom) panels. Odd-even, odd-odd 
and even-even concurrences approach at the factorizing field 
exactly the same side limits as those of figs. [2H3] respectively, 
which are here of the same order as the values outside this 
field. 



C ee are in fact maximum at the factorizing field. Note 
again that for b°/b e > 1, C ee > C oe in the vicinity of b°, 
despite the absence of even-even couplings. 

Finally, we depict in fig. [S] the site magnetizations 
M CTj = (s?) together with the concurrence C ai = Ci be- 
tween one spin and the rest of the chain (Eq. ([TU)) ). For 
a spin 1/2 chain with [H,P Z ] — 0, both quantities are 
strictly related, since the reduced density matrix for one 
spin in a state with definite parity is diagonal in the s z 



basis (pt 
hence 



2 sf 



as 



for p = x,y) and 



(30) 



Thus, Ci = 1 when (sf) = (zero field). At the fac- 
torizing field it approaches the side limits (fTU)) . The en- 
suing entanglement entropy can be evaluated as Si = 
-Y, v =±P»to&2Pv, with p± = |(1 ± y/l - Cf ). 

While for a uniform field the even and odd site mag- 
netizations coincide and decrease stepwise as the field 
increases, approaching —1/2 for strong fields, for non- 
uniform fields they first acquire opposite signs (M a = 
—M e ) in the "dimer phase", i.e., before the first par- 
ity transition. Here the magnetization is essentially de- 
scribed by the two-qubit result (|25|) . which yields M e = 
this symmetric system^. In the case considered C QO and — M a > in the state I^T) if b° > b e . Accordingly, 
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FIG. 5: (Color online) Magnetization at even and odd sites 
M a = (s%) (top panels) and the concurrence C a = Ci (Eq. 
(J30j) ) between the site and the rest of the chain (bottom pan- 
els), in the dimer chain of figs. [2}[3] with \ = 0.9. Results for 
an alternating field with fixed ratio b°/b e = 3 and for a uni- 
form field are depicted. The discontinuities at the factorizing 
field are explicitly shown. The "dimer phase" (fields below the 
first transition) presents opposite magnetizations and leads to 
a non-monotonous behavior of C e after this transition. 

M e first increases as b° (and hence b e and 6_ in (|25jl) 
increases, in close agreement with Eq. (|25j) . After the 
first transition, M e starts to decrease, crossing and ap- 
proaching — 1/2 (together with M a ) for strong fields, even 
though it may still increase between transitions. 

This entails a non-monotonous behavior of C e for in- 
creasing fields, particularly appreciable for low a, where 
C e saturates again (C e = 1) at a finite field, i.e., when 
M e vanishes. At the factorizing field C a approaches the 
limits (|21[) . which are independent of a, with the mag- 
netization step there given by^ AM; = (s?)_ — = 
bin ^(-ejei) ^ p Qr s ^ ron g figi^g the behavior of M a and 
C a can again be described perturbatively: We obtain 
M a M e » -|[1 - C-^j^n with C « C e « 

I 6 I" 



IV. CONCLUSIONS 



We have then examined the magnetic behavior of a 
finite spin 1/2 XY dimer chain. The factorizing field 
corresponds to the last parity transition exhibited by the 
exact GS for increasing field. Dimerization breakdown 
takes then place in steps, with all signatures of dimer 
structure being completely lost at the factorizing point: 
For a uniform field, the concurrence between any two 
spins approaches there (at each side) a constant value, 
independent of separation and coupling ratio v e /v°. The 
same behavior occurs in an alternating field, except that 
in this case there arc three different concurrence side lim- 
its at separability, which depend on the odd-even field ra- 
tio. The entanglement between spins unconnected by the 
coupling may here exceed that between connected pairs. 

The previous properties are not a particular feature of 
the system considered. For full range coupling, the same 
behavior is obtained at separability, as the eigenstates 
become there independent of the coupling range. The 
behavior of the concurrence between one spin and the rest 
of the chain has also been examined. An alternating field 
can induce opposite magnetizations at even and odd sites 
before the first transition, leading to a non-monotonous 
behavior of this concurrence for increasing fields, with 
two saturation points. The present results shed light on 
the complex behavior of entanglement in these systems, 
and its relation with factorization. The exposed features 
can make such finite critical systems of special interest 
for diverse applications. 

The authors acknowledge support from CIC (RR) and 
CONICET (NC, JMM) of Argentina. 



Appendix: Exact solution of the dimer chain with 
alternating field 

By means of the Jordan- Wigner transformation^, and 
for a fixed value p = ± of the global SVparity P z , we 
may exactly rewrite the dimer XY Hamiltonian (fT2")) as a 
quadratic form in standard fermion creation and annihi- 
lation operators c] , Cj . For an alternating field b 1 = b aj , 
with o~j — o, e the site parity, we obtain 



i-5 b°/v° 



We have first determined the factorization conditions 
for general dimer-type arrays with XYZ couplings in gen- 
eral transverse fields. We have also examined the entan- 
glement properties of the associated definite parity states, 
which constitute the actual GS side-limits at separability 
in a finite system, showing that weak but non-zero full 
range pairwise entanglement can be reached in the vicin- 
ity of factorizing fields. The possibility of an alternating 
and field dependent separable GS through a non-uniform 
field along a separability curve (Eq. (fTB")) ) has also been 
shown, for general spin. Border corrections to the field 
allow exact separability also in open chains. 



where n + 1 = 1, «± = j( v x ± v y) an d Vj = 1> Vj" = 
1 — 26 j n . By means of separate discrete parity dependent 
Fourier transforms for even and odd sites, 
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where k+ = {|, . . . , f — 5}, fc_ = {0 

12 12 x k( a L a ^ - 1) . 



v + — I 2 ' 

rewrite (jA.ll) as 



, § - 1}, we may and A = (b ° )2 + (b ° )2 + |^| 2 + \v k \ 2 . Care should be taken 
to select the correct signs of X k in order that the vac- 
, , , uum of the operators a kv has the proper SVparity and 

v + c ko c ke + v_c c _fc e + /i-c.J-^presents the lowest state for this parity. 



(A.2) 



kek„ v=± 



where v± — v± ±v±e lujk . The final diagonal form (|A.2|) 
is obtained by means of a Bogoliubov transformation 



ko 



Y,v=± u ka a ku + v ka- a ~ku determined through the 



diagonalization of 4 x 4 blocks 



n k = 



whose eigenvalues are ±AjJ", ±A fc , with 



/ 


b° 







-t£ \ 




—v'l 


b e 














-b° 


v\ 


V 


— 





r+ 


-b e ) 



(A.3) 



The spin correlations in the lowest states for each 
z-parity can then be obtained from the ensuing basic 
fermionic contractions fa — (c\cj) ~ 9ij = (c\cj), 
which can be directly obtained from the inverse Fourier 
transform of (c^c^,) = £ w V&i&, {c'{ a c'l ka ,) = 
Vko-U-ka'- We then obtain, through the use of 
Wick's theorem, (s 



— fiifjj f; 



and (sts T \ — 



^A 2 - \b°b e - {v\ + v k _){v k + - v k _)\ 2 (i-j)x(i- 



2. + „2 

= |[det(At.) ± det(Ar.)], where Af- are 



2 L" c n~ij/ V"y 

j) matrices of elements 2(/ 
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